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Abstract 

A selection of di-lepton events with significant missing transverse momentum has been 
performed using a total data sample of 237.4 pb~^ at e"'"e~ centre-of-mass energies of 
183 GeV and 189 GeV. The observed numbers of events — 78 at 183 GeV and 301 at 189 
GeV — are consistent with the numbers expected from Standard Model processes, which 
arise predominantly from W^W~ production with each W decaying leptonically. This 
topology is also an experimental signature for the pair production of new particles that 
decay to a charged lepton accompanied by one or more invisible particles. Discrimination 
techniques are described that optimise the sensitivity to particular new physics channels. 
No evidence for new phenomena is apparent and model independent limits are presented on 
the production cross-section times branching ratio squared for sleptons and for leptonically 
decaying charginos and charged Higgs. Assuming a 100% branching ratio for the decay 

— > i^Xi, where Xi is the lightest neutralino, we exclude at 95% CL: right-handed 
smuons with masses below 82.3 GeV for m~- — m^o > 3 GeV and right-handed staus with 

/i 

masses below 81.0 GeV for m-- — m^n > 8 GeV. Right-handed selectrons are excluded 

at 95% CL for masses below 87.1 GeV for m~- — m^o > 5 GeV, within the framework of 

e Ai 

the Minimal Supersymmetric Standard Model assuming /i < —100 GeV and tan/5 = 1.5. 
Charged Higgs bosons, H='=, are excluded at 95% CL for masses below 82.8 GeV, assuming 
a 100% branching ratio for the decay — > r^i/,-. 
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1 Introduction 



We report on a set of selected events containing two oppositely charged leptons and significant 
missing transverse momentum. Data are analysed from e"'"e~ collisions at LEP at centre-of-mass 
energies of 182.7 and 188.7 GeV with integrated luminosities corresponding to 56.4 pb~^ and 
181.0 pb~^, respectively. The number of observed events and their studied properties are found 
to be consistent with the expectations for Standard Model processes, which are dominated by 
the t^vl~V final state [l = e, /i, r) arising from W+W~ production in which both W bosons 
decay leptonically: W~ — >£~z7^. 

This topology is also an experimental signature for the pair production of new particles 
that decay to produce a charged lepton accompanied by one or more invisible particles, such 
as neutrinos or the hypothesised lightest stable supersymmetric [jl[ particle (LSP), which may 
be the lightest neutralino, Xi? or the gravitino, G. Experimentally, invisible particles may also 
be weakly interacting neutral particles with long lifetimes, which decay outside the detector 
volume. We present the results of searches for the following new particle decays: 

charged scalar leptons (sleptons): ^^Xi (o^' ~^ ^^G), where may be a selectron 

(e), smuon (/i) or stau (f) and is the corresponding charged lepton. 

charged Higgs: H"^ r^Ur- 

charginos: Xi i^u ("2-body" decays) or xf ^ i^ux^ ("3-body" decays). 

The search for charged scalar leptons provides constraints on the selectron mass and in- 
directly the electron-sneutrino mass (in models where these are related). These searches are 
therefore also relevant to interpreting the results of searches for chargino and neutralino pro- 
duction since the production cross-section and branching ratios depend on the slepton masses. 

In most respects the analysis is similar to our published searches at centre-of-mass energies 
of 161, 172 and 183 GeV 0, |^. The analysis is performed in two stages. The first stage consists 
of a general selection for all possible events containing a lepton pair plus missing transverse 
momentum (Section |^). In this context the Standard Model i'^i' i~T7 events are considered as 
signal in addition to the possible new physics sources. All Standard Model processes that do 
not lead to (.'^v i~V final states — e.g. e^e~i'^i~ and £^£"(7) — are considered as background 
and are reduced to a rather low level by the event selection. In the second stage the detailed 
properties of the events (e.g. the type of leptons observed and their momenta), which vary 
greatly depending on the type of new particles considered and on free parameters within the 
models, are used to separate as far as possible the events consistent with potential new physics 
sources from W"'"W~ and other Standard Model processes (Section ^). 

The other LEP collaborations have published searches for sleptons in this channel using 
data at <183 GeV §. 

In this paper we describe fully only those aspects in which the second stage of the analysis 
differs significantly from [Q]. These are: 

• The use of the acolinearity of the event as an additional likelihood variable, and use of the 
fact that the momentum distributions employed in the likelihood calculation vary with 
acolinearity. 
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• The use of an extended maximum likelihood technique to calculate the upper limits on 
the cross-section times branching ratio squared. 



2 OPAL Detector and Monte Carlo Simulation 

A detailed description of the OPAL detector can be found elsewhere 0. 

The central detector consists of a system of chambers providing charged particle tracking 
over 96% of the full solid angle inside a 0.435 T uniform magnetic field parallel to the beam axis. 
It consists of a two-layer silicon micro-strip vertex detector, a high precision drift chamber, a 
large volume jet chamber and a set of ^-chambers that measure the track coordinates along the 
beam direction. 

A lead-glass electromagnetic calorimeter is located outside the magnet coil. It provides, 
in combination with the forward detectors, which are lead-scintillator sandwich calorimeters 
and, at smaller angles, silicon tungsten calorimeters, geometrical acceptance with excellent 
hermeticity down to approximately 25 mrad. 

The magnet return yoke is instrumented for hadron calorimetry and consists of barrel and 
endcap sections along with pole tip detectors that together cover the region |cos^^| < 0.99. 
Outside the hadron calorimeter, four layers of muon chambers cover the polar angle range of 
|cos6'| < 0.98. Arrays of thin scintillating tiles have been installed in the endcap region to 
improve trigger performance, time resolution and hermeticity for experiments at LEP 2 
Of particular relevance to this analysis are the four layers of scintillating tiles (the MIP plug) 
installed at each end of the OPAL detector covering the angular range 43 < 6 < 220 mrad. 

The following Standard Model processes are simulated. 4-fermion production is simu- 
lated using the grc4f |^ generator at ^/s=183 GeV, and using the KoRALW [|| generator 
at ^/s=189 GeV. KoRALW uses the grc4f matrix elements to calculate the four-fermion cross- 
sections including interference effects and includes a detailed description of hard radiation from 
initial, intermediate and final state charged particles. Two-photon processes are generated using 



the program of Vermaseren [T^ and grc4f for e^e~£~^i~, and using Phojet [|ri[], Herwig [|T2 
and grc4f for e"'"e~qq. Because of the large total cross-section for e"'"e~e"'"e~, ji'^jj,^ and 
e"'"e~qq, soft cuts are applied at the generator level to preselect events that might possibly lead 
to background in the selection of t^vl'V final states. No generator level cuts are applied to 
the e^e~r"'"r~ generation. The production of lepton pairs is generated using Bhwide ||13[ and 



Teegg [III] for e+e"(7), and using KORALZ [ITH] for yU+/i~(7) and r+r'(7). The production of 



quark pairs, qq(g), is generated using Pythia and the final state vv^^ is generated using 
NuNUGPv98 IITTII and KoRALZ. 



Slepton pair production is generated using Susygen |T^. Chargino pair production is 



generated using Dfgt for three-body decays, and Susygen for two-body decays. Charged 



Higgs boson pair production is generated using HzHA and Pythia. 

All Standard Model and new physics Monte Carlo samples are processed with a full simula- 
tion of the OPAL detector and subjected to the same reconstruction and analysis programs 
as used for the OPAL data. 
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3 General Selection of Di-lepton Events with Significant 
Missing Momentum 



The general selection of acoplanar lepton pair events, which selects events containing low mul- 
tiplicity jets with significant missing transverse momentum, pf^^^^, is described in detail in 0. 
In we made use of the improved hermeticity for non-showering particles in the forward direc- 
tion provided by the MIP plug. Subsequent modifications have been made for the analysis of 
the data taken at 189 GeV, the most important of which was prompted by Monte Carlo studies 
which showed that, in the majority of Standard Model background events (non l^vl'V) ac- 
cepted by the general selection, p^^^^ was overestimated due to the mis-measurement of tracks 
and clusters. In the current analysis the uncertainty on p^^^^ is calculated from the measure- 
ment uncertainty on the observed tracks and clusters, and the requirement in ||^ that p™***^ 
exceed certain fixed cuts is replaced by the requirement that it exceed the cut values by at least 
one standard deviation of the calculated measurement uncertainty. These changes - both the 
MlP-plug cut and the p^^^^ significance, have markedly reduced the residual backgrounds and 
have allowed the selection efficiency to be substantially increased by removing many of the cuts 
of selection II which are now redundant. In brief, the essence of event selection II is now : 

• p™'^V-^beam should significantly exceed 0.05. 

• The scaled missing transverse momentum with respect to the transverse thrust axis, 
a™'^'^/i?beam should cxceed 0.022 for events with low acoplanarity. 

• Events with values of p™^^/-Ebeam which could potentially be balanced by beam energy 
muons in the MlP-plug acceptance are rejected if evidence for such forward-going particles 
is seen in these scintillators. 

The numbers of events passing the general selection at each centre-of-mass energy in the 
data are compared to the Standard Model Monte Carlo predictions in Table |I|. The total 
number of events predicted by the Standard Model is given, together with a breakdown into 
the contributions from individual processes. The number of observed candidates is consistent 
with the expectation from Standard Model sources, which is dominated by the t^v i~V final 
state arising mostly from W"*"W~ production in which both W's decay leptonically: W~ ^t~V(,. 



(GeV) 


data 


SM 






iiqq 


£+£-(7) 


z/i/77 


183 


78 


81.4±0.8 


77.5±0.7 


3.4±0.5 


0.07±0.03 


0.31±0.04 


0.06±0.03 


189 


301 


303.3±1.9 


292.5±1.6 


4.4±0.8 


1.3±0.1 


4.6±0.4 


0.46±0.04 



Table 1: Comparison between data and Monte Carlo of the number of events passing the 
general selection at each centre-of-mass energy. The total number of events predicted 
by the Standard Model is given, together with a breakdown into the contributions 
from individual processes. The Monte Carlo statistical errors are shown. 



The second stage of the analysis, in which we distinguish between Standard Model and 
new physics sources of lepton pair events with missing momentum, is described in Section ^. 
Discrimination is provided by information on the lepton identification, the acolinearity of the 
event, and the momentum and —qcos9 of the observed lepton candidates, where q and 9 are 



6 



Lepton 
identification 


V~s= 


189 GeV 




183 GeV 


data 


SM 


data 


SM 


e"'"e~ 


49 


45.2±0.7 


14 


12.1±0.3 




49 


47.0±0.7 


13 


13.7±0.3 




16 


11.0±0.5 


1 


2.5±0.2 




79 


83.4±0.9 


20 


25.6±0.4 


4- 1 -r 


26 


36.6±0.6 


8 


9.8±0.3 




40 


35.4±0.6 


8 


9.2±0.2 


e^, unidentified 


20 


19.6±0.7 


5 


3.8±0.2 


/i^, unidentified 


14 


17.8±0.5 


7 


3.7±0.2 


h^, unidentified 


8 


7.3±0.3 


2 


0.9±0.1 



Table 2: The lepton identification information in tlie events passing the general selection 
compared witli tie Standard Model Monte Carlo at each centre-of-mass energy. %" 
means that the lepton is identified neither as an electron nor muon and so is probably 
the product of a hadronic tau decay. Leptonic decays of taus are usually classified 
as electron or muon. ''Unidentified" means that only one isolated lepton has been 
positively identified in the event. 

the charge and polar production angle of the lepton. We check here on the degree to which 
these quantities are described by the Standard Model Monte Carlo. The lepton identification 
information in the event sample produced by the general selection at each centre-of-mass en- 
ergy is compared with the Standard Model Monte Carlo in Table |[ For the event sample at 
V^=189 GeV, Figure |I| shows the distributions of (a) the momentum scaled by the beam energy 
of each charged lepton candidate, (b) the value of —qcos9 of each charged lepton candidate, 
and (c) the acolinearity of the event. The data are compared with the Standard Model Monte 
Carlo predictions, which are dominated by the final state £+z/£~z7. 

The cuts used to veto two-photon background introduce an inefficiency in the event selec- 
tion due to random detector occupancy (principally in the SW, FD and MIP plug detectors) 
that is not modelled in the Monte Carlo. This inefficiency has been measured using randomly 
triggered events collected during normal data taking. For events with very low missing trans- 
verse momentum the inefficiency has a value of 8.3% at ^/s = 189 GeV and 8.2% at ^/s = 183 
GeV for events, and decreases to a negligible value for events with p™'^^/-Ebeam > 0.25. When 
quoting expected numbers of Standard Model events and selection efficiencies for potential new 
physics sources, the variation of veto inefficiency with p™'**** is taken into account. 

4 Likelihood Method to Distinguish Signal and Back- 
ground 

Starting from the general selection of acoplanar di-lepton events, we search for the production 
of new particles by using a likelihood technique which combines the information from a number 
of discriminating variables in order to distinguish between new physics signals and Standard 
Model sources of such events, the most important of which are i~^i'i~T7 and e'^e~£^i~. 

Discrimination between new particle pair production and the Standard Model background 
is performed by considering the likelihood that an event is consistent with being either signal 
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or background. Given an event, for which the values of a set of variables Xj are known, the 
likelihood, Ls, of the event being consistent with the signal hypothesis is calculated as the 
product of the probabilities, Ps{xi), that the signal hypothesis would produce an event with 
variable i having value Xi, Ls = Hi Ps{xi)- Similarly, the likelihood of an event being consistent 
with the background hypothesis is = IliPB{xi). The discriminating quantity used is the 
relative likelihood, L^, defined by: 



ijji = . 

Ls + Lb 

An event with high L/j is signal-like and an event with low L^ is background-like. 
The following quantities are used as likelihood variables (xj) in the analysis: 

(i) Scaled momentum, p/-Ebeam, of each lepton, 

(ii) Acolinearity of the event, defined as the supplement of the 3-D angle between the two 
leptons. 

(iii) —qcos9 for each lepton (smuons, staus and charged Higgs only), 

(iv) Lepton type variable (defined in Section |4.4|) . 

Sections [4.1| to [4.4| describe the properties of each of these variables for signal and back- 
ground. 



4.1 The Momentum Likelihood Variable 



The dominant Standard Model process leading to the acoplanar di-lepton signature arises from 
leptonically decaying W pairs, and the momentum distribution is highly populated between 
about 0.25 and 0.7 (Figure |l|(a)). The kinematics for the signal vary considerably with the 
mass difference. Am, between the parent particle (e.g. selectron) and the invisible daughter 
particle (e.g. Xi)^ since this determines how much energy is available to the lepton. The 
kinematics also vary to a lesser extent with the mass m of the parent particle, due to Lorentz 
boost effects. 

A significant change to the analysis with respect to earlier publications is the inclusion of 
the fact that for a given m and Am, the lepton momentum distribution varies according to 
the acohnearity of the event. Since the parent particles (e.g. sleptons) are produced back to 
back, then if an event has high acolinearity, one of the leptons will typically be travelling in a 
direction at an angle greater than 7r/2 to that of the parent particle (in the lab. frame). In 
this case, the lepton momentum in the lab. frame is reduced relative to its value in the rest 
frame of the parent particle and the lepton is therefore soft. An event with low acolinearity 
will in general have both leptons travelling in similar directions to the parent particles and the 
Lorentz boost results in both leptons having high momenta, provided that the parent particle 
mass is not close to the kinematic limit. Since the Lorentz boost is stronger for low parent 
particle mass, these effects are greater when m is small. 

Figure |^ shows momentum distributions in the {xmax,Xmin) plane, where Xmax and Xmin 
are the scaled momenta of the higher and lower momentum leptons, respectively, for a smuon 
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signal with m=45 GeV, Am=45 GeV, for 3 different ranges of acolinearity. Tlie corresponding 
plots for background are also shown. 

Signal and Standard Model Monte Carlos are used to construct reference histograms for 
momentum in a grid of points in the (m,Am) plane, with m ranging from 45 to 94 GeV, 
and Am ranging from 2 GeV to m. Each of the signal momentum distributions and the 
background distribution is subdivided into the following 3 ranges of acolinearity (in radians): 
< Oacoi < 0.8, 0.8 < 9acoi < 1-6 and 1.6 < 9acoi < tt. Each of these distributions is then further 
subdivided according to whether the observed lepton is the higher or lower momentum lepton. 

For the searches in which the final state particles can be the decay products of taus (staus, 
charginos and charged Higgs), each momentum probability distribution must be further sub- 
divided, as the momentum spectrum depends on the lepton identification. One momentum 
probability distribution is constructed for the case in which the observed lepton is identified 
as e or /i, and another for the case in which the observed lepton is identified as a hadronically 
decaying tau, or is unidentified. 

Probabilities Ps{xi) and PB{xi) for the likelihood calculation are found by reading from the 
appropriate reference histograms. 



4.2 The Acolinearity Likelihood Variable 

For signal, the distribution of the acolinearity angle varies with m and Am. For low parent 
particle mass, the Lorentz boost results in a tendency for the leptons to be in the directions 
of the parent particles, resulting in the acolinearity being peaked towards low values, whereas 
for high mass, the parent particles are produced close to being at rest, and the leptons have no 
preferred direction. For background, the acolinearity distribution is peaked towards low values 
due to the spin structure of the weak couplings. 

Figure |^ shows some example acolinearity distributions for signal, which can be compared 
to the distributions for background and data shown in Figure |I|. 

The use of the acolinearity as a likelihood variable is complementary to its use in defining 
the momentum probability, in that the masses for which it offers the greatest discrimination 
as a likelihood variable are the masses where the gain in distinguishing power described in 
Section |4.1| is small, and vice versa. 



4.3 The —qcos9 Likelihood Variable 

As described in f^, the distribution of the quantity —qcos9 , where q and 9 are the charge and 
production angle of an observed lepton, is forward peaked for W"'"W~ production due to the 
dominance of the neutrino exchange amplitude and the V-A nature of W decay, whereas for 
smuon, stau and charged Higgs production the distribution is symmetric and peaked towards 
I cos6'| = 0, due to the scalar nature of these particles. 

This variable is not used in the likelihood calculation for selectrons or charginos because 
these particles can be produced via t-channel neutralino exchange and sneutrino exchange, 
respectively, in addition to s-channel production. This results in the expected —qcos9 distri- 
bution of selectrons and charginos being model-dependent and potentially similar to that of 
the W^W~ background. 
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4.4 The Lepton Type Likelihood Variable 



A value is assigned to the lepton type variable according to which types of lepton are identified 
in the event. There are nine possible values, corresponding to the nine event types listed in 
table |. 

For the selectron (smuon) analysis a cut is applied at the same time as the general selection, 
requiring at least one electron (muon) and no muons (electrons), reducing the background by 
about 2/3 (depending on slepton type and a/s) with negligible loss of efficiency. With this cut 
applied, there are only three possible values of the lepton type variable. 



4.5 Lji and Lb Distributions 

For each search channel, reference histograms are constructed for each of the likelihood variables 
at each point in m and Am for which signal Monte Carlo has been generated. A smoothing 



algorithm [|22[ is applied to the histograms to reduce the effects of statistical fluctuations. The 
reference histograms are then used to construct La distributions. 

distributions for signal Monte Carlo, Standard Model Monte Carlo and data are shown 
in Figure ^ for the speciflc example of the analysis for smuons with a mass of 80 GeV and a 
smuon-neutralino mass difference of 60 GeV. There is considerable variation in the shapes of 
these distributions with m and Am. 

A check of consistency between data and the Standard Model can be performed without 
reference to a particular signal by comparing the L b distributions for data and Standard Model. 
Figure |(a) shows the Lb distributions for the Standard Model (histogram) and data (points) 
for events passing the general selection. All the likelihood variables are used. Figures ^(b) 
and (c) show the same information after making the initial lepton identiflcation requirements 
given in Section [4.4| for the selectron and smuon searches respectively. In Figure ^(b), only 



the variables used in the selectron analysis are used. In each of the plots, the secondary peak 
at high Lb corresponds to events which have only one identifled lepton and therefore fewer 
variables entering the likelihood [[ In all three plots, the data is in good agreement with the 
Standard Model. 



5 Calculation of Cross-Section Limits 
5.1 Introduction 

In 1^ , the limit on the cross-section was calculated by flnding an optimised cut on the value of 
Lji as a function of m and Am for each centre of mass energy, and applying this cut to signal 
Monte Carlo, Standard Model Monte Carlo and data. The resulting efficiencies, expected 
backgrounds and numbers of candidates were used to calculate cross-section limits using the 



likelihood ratio method [23 to combine the information. 



In this paper, we describe the use of an extended maximum likelihood calculation to deter- 
mine the cross-section limits. In this method, no cut is applied on L/j. Information contained 



^This effect cancels when the likehhood ratio L^j is calculated because the events will have high Ls for the 
same reason. 
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in the values of each individual candidate event, and in the shapes of the Lr distributions 
for signal and background are used as input to the limit calculation, rather than the numbers 
of events passing a cut. In this way, considerably more of the available information is used. 

The advantage of a cut free method can be seen by considering, for example, a case where 
there is an excess of candidates passing a cut on L^,- The information about whether the events 
all lie close to the cut, or whether they are clustered towards Lr=\ (suggesting the presence of 
a signal) is not used. The use of the additional information makes the analysis more sensitive 
for discovery, and at the same time is able to set more stringent limits in the absence of signal. 
The expected sensitivity (ie., the expected upper limit on the cross-section) is improved by as 
much as 20%, depending on m and Am, using this technique. 

5.2 Extended Maximum Likelihood Technique 

The upper limit on the cross-section times branching ratio squared at 95% confidence level, 
(795, is calculated by forming a likelihood, iv((7s), of the set of Lr values for the data being 
consistent with the expected Lr distribution for Standard Model plus a signal produced with 
cross-section times branching ratio squared as- ergs is the value of cr^ below which 95% of the 
area under the likelihood function lies. 

5.2.1 The Likelihood Function 

Extended maximum likelihood combines standard maximum likelihood with the Poisson prob- 
ability of observing N candidate events when v are expected: 

-V N N 
■ i=l 

where P{Lji.; B, S) is the probability of event i having Lr = Lr^, given Lr distributions B 
and S for background and signal. 

Dropping the constant N\, this can be re- written: 

N 

lnL = -iy + J2HQiLR,;B,S)] 
1=1 

where Q is identical to P but normalised to u instead of 1 (Q = I'P). 
The expected number of candidates u is given by: 

where /is is the expected number of Standard Model events with non-zero Lr passing the 
general selection (similarly, N is the number of data candidates with Lr 7^ 0), e is the signal 
selection efficiency of the general selection, C is the experimental luminosity and a; is a weight 
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factor which takes into account that the expected production cross-section varies with ^/s, but 
the hmit on the observed cross-section is quoted at -y/s = 189 GeV. 



where (Jigg is the expected cross-section for ^/s = 189 GeV and CTj is the expected cross-section 
for the i'th value of y^. For scalar particles, for example sleptons, we assume that the expected 
cross-section varies as For spin | particles, for example charginos, we assume that the 

expected cross-section varies as j3/ s. 

The function Q is the probability of event i having Lji = L^. , given Lr distributions B and 
S for background and signal, normalised to u. This is given by: 

Q = ij,bB{Lr^) + eCLLiasS{LR^), 

where the functions B and S, formed using background and signal Monte Carlo respectively, 
are normalised to 1. 

Hence the likelihood function is given by: 

N 

lnL{as) = -{jJiB + eCuGs) ^ln[/XB5(Lij.) -|- eCuasS{Lii.)]. 

i=l 



5.2.2 Limit Calculation 

For data at a single centre-of-mass energy, the upper limit on the cross-section at 95% confidence 
level is the value of (795 which satisfies: 



/„~L(<7,)d<T. 

The generalisation to Necm values of ^/s is: 

189 

where cr^^^ is the cross-section at \/s = 189 GeV. 



5.2.3 Verification of the Method 

The technique used to calculate limits was tested using a toy Monte Carlo to simulate data 
sets for an ensemble of simulated experiments in which a signal is present with cross-section 

as- The total number of candidates was drawn from a Poisson distribution with mean v and 
for each candidate, a value of Lr was assigned, chosen randomly according to the sum of the 
expected Lr distributions for background and signal with cross-section Ug- 
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(jgs was calculated for 500 simulated experiments. This was done using smuon signals at all 
m and Am for which Monte Carlo has been generated, for 189 GeV separately and for 183 GeV 
and 189 GeV combined, and for a number of values of ag- In all cases, 95% of the 500 (T95 values 
were found to be greater than the true cross-section as, to within the statistical error expected 
from the finite number of simulated experiments. 

As a further test, the best estimate of the signal cross-section, abest, was calculated for each 
simulated experiment. This is the value of as at which the likelihood function Li{as) peaks. 
The abest distribution from the 500 simulated experiments was in all cases found to peak at the 
true cross-section. 

5.3 Limit Calculation at an Arbitrary Point in m and Am 

Monte Carlo signal events are available only at certain particular values of m and Am. The 
values of m range typically from m = 45 GeV up to m «i i?beam in 5 GeV steps. The values of 
Am vary between 2 and Am= m. In order to calculate (Xgs at an intermediate point in m and 
Am, it is necessary to be able to calculate at that point for a given event, which requires 
the existence of reference histograms for the likelihood variables for any m and Am. 

An algorithm has been developed to construct the reference histograms at any intermediate 
value of m and Am, given the histograms at the four nearest signal Monte Carlo grid points, 
assuming a linear variation in the shape of the histograms with m and Am. This procedure 
has been tested by re-constructing histograms at gridpoints using the histograms at adjacent 
gridpoints. 

For a given point m and Am, the signal and background distributions and the data Ln 
values are calculated using the interpolated reference histograms. Signal Monte Carlo at an 
intermediate point is simulated using signal Monte Carlo events at the nearest grid point. This 
is done by re-defining the value of each likelihood variable for an event, such that the fraction of 
the corresponding reference histogram at the intermediate point which lies below the re-defined 
value is the same as the fraction of the histogram at the grid point which lies below the original 
value. 

The effects of statistical fluctuations in the distributions are reduced using the same 
smoothing algorithm as applied to the reference histograms (section [4.5|) . 

The remaining input to the limit calculation is the signal efficiency, obtained at intermediate 
values of m and Am by linear 2-dimensional interpolation. 

6 New Particle Search Results 

We present limits on the pair production of charged scalar leptons, leptonically decaying charged 
Higgs bosons and charginos that decay to produce a charged lepton and invisible particles. 

The 95% CL upper limit on new particle production at y/s = 189 GeV, obtained by combin- 
ing the data at y/s = 189 GeV and = 183 GeV is calculated at each kinematically allowed 
point on a 0.5 GeV by 0.5 GeV grid of m and Am, using the Lji distributions for signal and 

^Particle masses less than 45 GeV are not considered because these masses were accessible at LEPl and 
because radiative return to the Z means that the event topology can be different. 
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background, the values of the data events, and the efficiency of the general selection at that 
point as input. 

In addition to the Monte Carlo statistical error on the signal efficiency, we assign a 10% 
systematic error on the estimated selection efficiency to take into account uncertainties in 
trigger efficiency, detector occupancy, lepton identification efficiency, luminosity measurement, 
the interpolation procedure, and deficiencies in the Monte Carlo generators and the detector 
simulation. An additional systematic error in the stau analysis is the effect of tau polarisation 
in the modelling of the stau signal. It is possible for the tau produced in stau decay to have 
any polarisation value in the range [-1,1] |2^. This was studied by using stau Monte Carlo 



events with tau polarisations of +1, and —1 to determine the amount by which the expected 
limit on the cross-section times branching ratio squared is overestimated or underestimated if 
a polarisation of zero is assumed when the true polarisation is +1 or —1. The size of this effect 
was found to vary with m and Am, but to be always less than about 5% , and so is included 
in the 10% systematic error. 

At high values of Am the dominant background to the searches for new physics results from 
W"'"W~ production. High statistics Monte Carlo samples for this process are available that 
describe well the OPAL data [Q. In addition to the Monte Carlo statistical error, we assign 
a 10% systematic error on the estimated background to take into account uncertainties in the 
shapes of the L/j distributions and reference histograms, and in the interpolation procedure, 
and deficiencies in the Monte Carlo detector simulation. At low values of Am the dominant 
background results from e'^e~i'^i~ events. The background uncertainty at low Am is dominated 
by the limited Monte Carlo statistics; the uncertainty is typically 20-80% at low Am. In setting 
limits the Monte Carlo statistical errors and other systematics are taken into account according 



to the method described in 26 



6.1 Limits on Production Cross-section Times Branching Ratio Squared 

Limits on the production cross-section times branching ratio squared for new physics processes 
are now presented in a manner intended to minimise the number of model assumptions. The 
95% CL upper limits at y/s = 189 GeV shown in Figures ^ - |ll] are obtained by combining 
the data at the two centre-of-mass energies 183 and 189 GeV using the assumption that the 
cross-section varies as for sleptons and f3/s for charginos. The chosen functional forms 
are used for simplicity in presenting the data and represent an approximation, particularly for 
processes in which t-channel exchange may be important, that is, selectron pair and chargino 
pair production. In these cases the cross-section dependence on centre-of-mass energy is model 
dependent, depending on the mass of the exchanged particles and the couplings of the neutrali- 
nos and charginos. The selectron Monte Carlo events were generated at /i = —200 GeV and 
tan/5 = 1.5 using Susygen. We have found by varying fi and tan/3 that the above choice gives 
a conservative estimate of the selection efficiency for selectrons. 

Upper limits at 95% CL on the selectron pair cross-section at ^/s = 189 GeV times branch- 
ing ratio squared for the decay e~— ^ e~Xi ^ire shown in Figure ^ as a function of selectron 
mass and lightest neutralino mass. These limits are applicable to SlGl and e^e^ production. 
The corresponding plots for the smuon and stau pair searches are shown in Figures ^ and |^, 
respectively. Note that if the LSP is the gravitino, G (effectively massless), then for prompt 
slepton decays to a lepton and a gravitino the experimental signature would be the same as 
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that for i — i>£ with a massless neutrahno. In this case the hmits given in Figures |^ - ^ for 
TTL^o = may be interpreted as hmits on the decay i~ ^ i~G 

The upper hmit at 95% CL on the chargino pair production cross-section times branching 
ratio squared for the decay xf £^ue (2-body decay) is shown in Figure |^. The hmit has been 
calculated for the case where the three sneutrino generations are mass degenerate. The upper 
limit at 95% CL on the chargino pair production cross-section times branching ratio squared 
for the decay Xi^ W^Xi ~^ ^^i^Xi (3-body decay) is shown in Figure |T^. 

The upper limit at 95% CL on the charged Higgs boson pair production cross-section times 
branching ratio squared for the decay H"^ r^z/,- is shown as a function of mH+ as the solid line 
in Figure |lT[ The limit is obtained by combining the 183 and 189 GeV data-sets assuming the 
mH+ and ^/s dependence of the cross-section predicted by Pythia, which takes into account 
the effect of initial state radiation. The dashed line in Figure [Tl| shows the prediction from 
Pythia at ^/s = 189 GeV for a 100% branching ratio for the decay H='= — >• t^Ut- With this 
assumption we set a lower limit at 95% CL on mH+ of 82.8 GeV. 

6.2 Expected Limits and Confidence Levels for Consistency with 
Expectation 

Table |] gives the values of the following quantities for a number of values of m and Am in the 
search for selectrons: 

1. The signal selection efficiency of the general selection at 189 GeV (the efficiencies at 
183 GeV are similar). 

2. The 95% CL upper limit on the cross-section times branching ratio squared at 189 GeV, 
obtained by combining the data at ^/s = 189 GeV and ^/s = 183 GeV. 

3. The expected 95% CL upper limit on the cross-section times branching ratio squared 
in the absence of signal (cxgs). This is calculated using an ensemble of 1000 toy Monte 
Carlo experiments to simulate the data, in which the total number of candidates for each 
experiment is drawn from a Poisson distribution with mean equal to the number of events 
expected from the Standard Model, and for each candidate, a value of is assigned, 
chosen randomly according to the expected Ln distribution for Standard Model processes. 
The expected limit at a given point in m and Am is the mean value of the limit for the 
ensemble of simulated experiments. 

4. The confidence level for consistency with the Standard Model, calculated as the fraction of 
the simulated experiments for which the upper limit on the cross-section times branching 
ratio squared is greater than or equal to the value calculated using the actual data. In 
the absence of signal, a CL of 50% is expected on average 0. 

Tables H, ^, H, ^ and H show the same information for smuons, staus, charginos with two-body 
decay, charginos with three-body decay and charged Higgs, respectively. 

^Values of 100% correspond to points where there are no candidate events with non-zero in the OPAL 
data. In this case, all toy Monte Carlo experiments will have a value of (T95 equal to or (if there are candidates 
with non-zero Lji) greater than the value for the actual data. 
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For some points in m and Am in Tables ^ and ^, the confidence level for consistency with 
the Standard Model is small (around 1%). The probability of getting a low confidence level 
for one or more points in m and Am for one or more of the search channels depends on the 
degree of correlation among the different (m, Am) points and among the different channels. 
The degree of correlations between adjacent points is strong when the momentum distributions 
for those points are similar. The momentum distributions vary slowly with both m and Am 
when Am is high (hence the clustering of low confidence level values in Table but vary 
considerably with Am when Am is low. 

This effect was investigated by calculating the cross-section limits for each of 1000 experi- 
ments in which the data is simulated by randomly selected Standard Model Monte Carlo events. 
For each experiment, the number of events taken from a Monte Carlo sample simulating a given 
process is drawn from a Poisson distribution with mean equal to the number of events expected 
for that process. For each experiment, the confidence level at each (m. Am) point at which 
signal Monte Carlo has been generated was calculated as already described, and the number of 
experiments for which a confidence level of 0.6% Q or less is obtained for at least one point in 
m and Am in at least one search channel was determined. This was found to be the case for 
390 of the 1000 experiments. 

As a cross-check, taking the mean of the ensemble of limits obtained at each (m. Am) point 
for these simulated experiments was used as an alternative to the method described above to 
obtain ((T95). The results were found to be consistent. 



6.3 Limits on New Particle Masses 

We can use our data to set limits on the masses of right-handed sleptonsQ based on the expected 
right-handed slepton pair cross-sections and branching ratios. The cross-sections have been 
calculated using Susygen at each centre-of-mass energy and take into account initial state 
radiation. In Figure ^ we show limits on right-handed smuons as a function of smuon mass 
and lightest neutralino mass for several assumed values of the branching ratio squared for 
/i^ fi^Xi- The expected limit, calculated using Monte Carlo only, for a branching ratio of 
100% is also shown. For a branching ratio fl^ /^^Xi of 100% and for a smuon-neutralino 
mass difference exceeding 3 GeV, right-handed smuons are excluded at 95% CL for masses 
below 82.3 GeV. The 95% CL upper limit on the production of right-handed f+r^ times the 
branching ratio squared for t^Xi is shown in Figure |l^. The expected limit for a 

branching ratio of 100% is also shown. For a branching ratio '^^Xi of 100% and for a 

stau-neutralino mass difference exceeding 8 GeV, right-handed staus are excluded at 95% CL 
for masses below 81.0 GeV. No mixing between and tr is assumed. However, the cross- 
section ratio a-+-- /a-+-^ at ^/s=189 GeV varies between 0.89 and 1.20, depending only on 
the mixing angle. Using this information, the limits shown in Figure |l3| can be applied to 
any degree of stau mixing by multiplying the predicted cross-section for f^f^ by the value of 



a-+--/a-+-- corresponding to the mixing angle considered. The hatched region in Figure |T3 

11 R R 

shows the range of possible positions of the line defining the excluded region for a branching 
ratio ff" — >• r^Xi of 100% for any degree of stau mixing. 

^This is the lowest value of the confidence level in Tables || to ^. 

^ The right-handed slepton is expected to be lighter than the left-handed slepton. The right-handed one 
tends (not generally valid for selectrons) to have a lower pair production cross-section, and so conventionally 
limits are given for this (usually) conservative case. 
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For the case of a massless neutralino (or gravitino) and 100% branching ratio, right-handed 
smuons and staus are excluded at 95% CL for masses below 85.4 GeV and 81.1 GeV, respec- 
tively, and ff" is excluded at 95% CL for masses below 80.0 GeV, for any degree of stau mixing. 

An alternative approach is to set limits taking into account the predicted cross-section 
and branching ratio for specific choices of the parameters within the Minimal Supersymmetric 
Standard Model (MSSM)^. For /i < -100 GeV and for two values of tan/? (1.5 and 35), 



Figures |Tj, |T5| and |T6| show 95% CL exclusion regions in the {'mi±, m^o) plane for right- 
handed selectrons, smuons and staus, respectively. For fi < —100 GeV and tan/3 = 1.5, right- 
handed sleptons are excluded at 95% CL as follows: selectrons with masses below 87.1 GeV for 

m~- — m^o > 5 GeV; smuons with masses below 81.7 GeV for m~~ — m^o > 3 GeV: and staus 
e ^1 /I ^1 

with masses below 75.9 GeV for m~- — m^o > 7 GeV. 



7 Summary and Conclusions 

A selection of di-lepton events with significant missing transverse momentum is performed 
using a total data sample of 237.4 pb~^ at e"'"e~ centre-of-mass energies of 183 and 189 GeV. 
The observed numbers of events, 78 at 183 GeV and 301 at 189 GeV, are consistent with the 
numbers expected from Standard Model processes, dominantly arising from W"'"W~ production 
with each W decaying leptonically. 

Further discrimination techniques are employed to search for the pair production of charged 
scalar leptons, leptonically decaying charged Higgs bosons and charginos that decay to produce 
a charged lepton and invisible particles. No evidence for new phenomena is apparent and model 
independent limits on the production cross-section times branching ratio squared for each new 
physics process are presented. 

Assuming a 100% branching ratio for the decay £^ — > i'^Xi, we exclude at 95% CL: right- 
handed smuons with masses below 82.3 GeV for — iti^q > 3 GeV and right-handed staus 
with masses below 81.0 GeV for m~_ — m^o > 8 GeV. Right-handed selectrons are excluded 

at 95% CL for masses below 87.1 GeV for m~_ — m^o > 5 GeV within the framework of the 

e Ai 

MSSM assuming < —100 GeV and tan/5 = 1.5. Charged Higgs bosons are excluded at 95% 
CL for masses below 82.8 GeV, assuming a 100% branching ratio for the decay —>■ t'^Vt-. 

The cross-section times branching ratio squared limits from the selectron, smuon and two- 
body chargino searches presented here are used in the interpretation of the results of |^ in 
terms of mass limits on charginos and neutralinos. 
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Am 








;GeV) 






(GeV) 


45 


55 


65 


75 


85 


94 


signal selection efficiency of the 


general selection at 189 GeV (%) 




2 


13. Oil. 1 


10.5±1.0 


7.8±0.8 


4.2±0.6 


1.4±0.4 


0.4±0.2 




25.1±1.4 


21.4±1.3 


22.1±1.3 


18.5±1.2 


12.9±1.1 


7.4±0.8 


K 
o 


tJLI.U— 1— -L •\J 


59.5±1.6 


59.0±1.6 


60.5±1.5 


en o_|_i K 


55.9±1.6 


1 n 


71. Oil. 4 


74.4±1.4 


76.7±1.3 


76.3±1.3 


I I .O -i — L . O 


76.0±1.4 


20 


78.8±1.3 


81.6±1.2 


85.2±1.1 


85.2±1.1 


84.6±1.1 


85.9±1.1 


m/2 


78.5±1.3 


84.8±1.1 


87.4±1.0 


90.6±0.9 


90.4±0.9 


92.1±0.9 


m- 20 


79. Oil. 3 


85.6±1.1 


88.8±1.0 


90.4±0.9 


91.8±0.9 


93.1±0.8 


m- 10 


78.9±1.3 


84.9±1.1 


89.3±1.0 


90.3±0.9 


91.1±0.9 


92.5±0.8 


m 


77.4±1.3 


84.1±1.2 


88.7±1.0 


90.1±0.9 


91.3±0.9 


93.2±0.8 


95% CL 


upper limit 


on the cross-section times BR^{c - 


-cx\0 (fb) 




2 


93.8 


117.2 


164.8 


305.5 


973.8 


4138.6 


2.5 


49.2 


56.5 


56.1 


67.7 


106.1 


223.7 


5 


36.7 


25.7 


24.3 


21.9 


23.3 


29.6 


10 


63.1 


50.2 


43.8 


30.8 


26.8 


22.0 


20 


64.8 


50.9 


44.7 


48.9 


54.7 


25.6 


m/2 


92.3 


73.2 


79.4 


54.5 


35.3 


20.0 


m- 20 


97.8 


101.5 


93.3 


87.8 


63.4 


23.5 


m- 10 


91.5 


76.9 


73.8 


89.9 


70.4 


29.6 


m 


81.4 


84.6 


85.8 


86.8 


64.0 


34.5 


expected upper limit on the cross-section times _Bi?^ (e 


- ex^) (fb) 


2 


148.3 


175.1 


237.0 


414.5 


1325.7 


4616.9 


2.5 


82.0 


90.9 


87.4 


104.9 


136.6 


264.2 


5 


42.2 


39.3 


38.9 


34.6 


33.8 


34.1 


10 


40.4 


33.3 


28.5 


25.8 


24.5 


24.8 


20 


69.1 


52.8 


43.2 


34.8 


27.9 


23.6 


m/2 


76.1 


71.5 


62.5 


62.0 


57.2 


33.6 


m- 20 


80.7 


78.7 


79.3 


80.0 


70.3 


38.9 


m- 10 


84.9 


78.6 


79.6 


81.1 


69.3 


37.9 


m 


81.9 


76.0 


76.0 


80.8 


71.1 


38.8 


CL for consistency with SM (%) 


2 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


2.5 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


5 


57.0 


88.6 


90.5 


95.4 


89.1 


100.0 


10 


5.2 


6.8 


6.3 


18.0 


33.1 


83.9 


20 


47.6 


44.6 


38.2 


10.4 


0.6 


19.7 


m/2 


20.7 


37.7 


17.9 


53.6 


88.8 


99.9 


m- 20 


20.3 


14.4 


22.8 


27.8 


52.0 


95.0 


m- 10 


30.5 


43.2 


45.6 


27.4 


37.2 


72.3 


m 


38.2 


28.1 


24.9 


30.9 


51.2 


56.0 



Table 3: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 

on the cross-section times BR^{e — > ex?), expected upper limit on the cross-section times 
BB?{e exi), and the confidence level for consistency with the Standard Model in the search 
for e+e" production for different values of m~- and Am. 



20 



Am 






m~- 


(GeV) 






(GeV) 


45 


55 


65 


75 


85 


94 


signal selection efficiency of the 


general selection at 189 GeV (%) 




2 


14.9±1.1 


14.1±1.1 


9.8±0.9 


6.5±0.8 


0.7±0.3 


O.OiO.O 


2.5 


26.7±1.4 


27.7±1.4 


25.0±1.4 


21.3±1.3 


14.3±1.1 


8.0±0.9 


5 


58.7±1.6 


60.5±1.5 


60.4±1.5 


60.0±1.5 


60.3±1.5 


57.4±1.6 


10 


75.8±1.4 


76.7±1.3 


76.2±1.3 


76.4±1.3 


76.5±1.3 


74.1±1.4 


20 


85.4±1.1 


86.0±1.1 


84.8±1.1 


84.2±1.2 


84.0±1.2 


83.2±1.2 


m/2 


86.3±1.1 


87.8±1.0 


87.2±1.1 


90.4±0.9 


90.4±0.9 


91.7±0.9 


m- 20 


86.7±1.1 


89.2±1.0 


89.7±1.0 


91.2±0.9 


92.9±0.8 


92.8±0.8 


m- 10 


88.5±1.0 


90.6±0.9 


89.6±1.0 


91.3±0.9 


92.4±0.8 


93.0±0.8 


m 


88.6±1.0 


90.6±0.9 


89.8±1.0 


91.1±0.9 


92.1±0.9 


92.5±0.8 


95% CL 


upper limit 


on the cross-section times Bli^^fl 


-MX^) (fb) 


2 


105.6 


119.3 


177.3 


300.8 


2069.3 




2.5 


63.3 


61.5 


72.0 


88.8 


98.8 


206.9 


5 


30.3 


31.0 


39.5 


31.8 


34.1 


28.8 


10 


35.0 


22.7 


18.5 


17.3 


18.2 


22.3 


20 


53.5 


38.7 


39.9 


41.3 


37.6 


25.5 


m/2 


64.7 


76.7 


81.6 


71.1 


51.7 


44.1 


m- 20 


80.1 


100.5 


104.6 


89.6 


45.3 


34.3 


m- 10 


85.4 


94.6 


92.8 


62.9 


43.4 


34.3 


m 


83.5 


87.7 


88.8 


58.9 


43.5 


34.7 


expected upper limit on the cross-section times {fl 


^ A^X^) (fb) 


2 


112.8 


116.1 


169.0 


286.8 


2977.5 




2.5 


67.6 


64.7 


68.9 


82.4 


119.5 


226.1 


5 


34.8 


33.8 


31.2 


28.2 


28.8 


30.2 


10 


30.8 


27.4 


25.3 


22.6 


22.6 


24.3 


20 


47.5 


39.3 


35.9 


31.3 


25.3 


22.1 


m/2 


50.9 


50.9 


50.4 


51.0 


50.1 


33.4 


m- 20 


51.8 


52.8 


55.8 


58.7 


56.8 


36.4 


m- 10 


53.1 


51.1 


53.7 


56.5 


56.5 


38.1 


m 


50.6 


49.7 


52.9 


57.6 


56.7 


39.5 


CL for consistency with SM (%) 


2 


64.6 


47.9 


45.7 


43.7 


100.0 




2.5 


65.7 


60.6 


54.1 


40.2 


100.0 


100.0 


5 


57.8 


47.5 


19.0 


35.9 


27.2 


100.0 


10 


26.1 


63.8 


76.1 


100.0 


100.0 


100.0 


20 


27.8 


42.1 


30.0 


15.0 


7.5 


8.5 


m/2 


17.6 


7.6 


4.8 


10.6 


36.0 


14.2 


m- 20 


6.3 


1.7 


1.1 


6.1 


66.8 


49.0 


m — 10 


4.4 


1.3 


2.7 


29.1 


69.0 


52.3 


m 


3.4 


2.8 


3.9 


36.3 


71.3 


57.1 



Table 4: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 

on the cross-section times BR^{fi fJ'Xi): expected upper limit on the cross-section times 
BR^{fl /iXi)) the confidence level for consistency with the Standard Model in the search 
for fl~^fl~ production for different values of m~- and Am. 



21 



Am 






m~- 


(GeV) 






(GeV) 


45 


55 


65 


75 


85 


94 


signal selection efficiency of the 


general selection at 189 GeV (%) 




2 




O.libO.O 


O.libO.O 


O.OibO.O 


O.OibO.O 


O.OibO.O 




7+0 1 


0.5±0.1 


0.3±0.0 


0.3±0.0 


O.libO.O 


O.OibO.O 


K 
o 


J. (J . U — 1— W . tJ 


14.0±0.3 


13.1±0.3 


11.6±0.3 


10.3ib0.3 


9.1ib0.3 


1 n 


QO o_|_n 4 


39.1±0.4 


39.5±0.4 


38.8±0.4 


39.4ib0.4 


38.8ib0.4 


20 


57.4±0.4 


59.4±0.4 


59.9±0.4 


60.9±0.4 


60.9±0.4 


62.2ib0.4 


m/2 


59.3±0.4 


64.8±0.4 


69.6±0.4 


71.8±0.4 


73.7ib0.4 


74.5ib0.4 


m- 20 


60.6±0.4 


69.1±0.4 


73.4±0.4 


74.4±0.4 


77.5ib0.4 


78.3ib0.4 


m- 10 


65.2±0.4 


71.1±0.4 


73.9±0.4 


76.2±0.4 


77.7ib0.4 


79.2ib0.4 


m 


66.2±0.4 


71.3±0.4 


74.5±0.4 


76.1±0.4 


77.7ib0.4 


79.0ib0.4 


95% CL 


upper limit 


on the cross-section times BR^{f 


- rxl) (fb) 


2 


12180.7 


32560.6 


16554.6 








2.5 


2678.7 


4145.1 


5526.8 


6808.6 


19454.0 




5 


137.6 


168.4 


181.7 


243.9 


274.5 


324.1 


10 


106.8 


92.1 


81.5 


75.7 


74.2 


86.7 


20 


101.1 


92.8 


84.7 


86.5 


77.1 


59.2 


m/2 


100.8 


97.5 


87.4 


76.3 


76.2 


90.6 


m- 20 


106.4 


93.5 


78.7 


75.8 


76.0 


87.6 


m- 10 


119.8 


97.1 


83.2 


75.5 


73.7 


97.6 


m 


115.4 


94.1 


86.1 


76.6 


81.3 


87.6 


expected upper limit on the cross-section times BR^{f 


-^rx\) (fb) 


2 


11163.9 


20013.4 


19897.3 








2.5 


2916.5 


4380.8 


6850.2 


7709.0 


17774.8 




5 


189.5 


209.3 


221.8 


237.2 


262.6 


284.6 


10 


107.0 


103.4 


97.4 


98.2 


94.2 


92.1 


20 


103.4 


95.8 


89.7 


82.9 


79.1 


74.3 


m/2 


106.9 


104.7 


99.5 


95.1 


93.8 


96.1 


m- 20 


108.3 


107.1 


108.3 


107.5 


108.1 


112.2 


m- 10 


113.9 


112.1 


114.4 


109.7 


114.4 


111.6 


m 


118.2 


114.1 


112.4 


110.6 


113.7 


114.0 


CL for consistency with SM (%) 


2 


37.4 


7.4 


100.0 








2.5 


56.8 


41.2 


77.2 


65.5 


45.3 




5 


75.5 


66.6 


62.8 


38.6 


32.8 


23.0 


10 


40.7 


56.8 


61.3 


70.0 


66.6 


50.3 


20 


43.5 


44.0 


46.9 


37.7 


42.5 


67.2 


m/2 


46.1 


47.5 


54.7 


64.4 


64.1 


47.5 


m- 20 


41.4 


54.3 


74.1 


76.0 


75.0 


67.0 


m- 10 


33.0 


56.0 


72.9 


79.0 


83.2 


54.0 


m 


40.6 


59.7 


67.4 


79.0 


73.3 


69.6 



Table 5: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 

on the cross-section times BF?{f — >■ rx?), expected upper limit on the cross-section times 
BR^ij TXi), and the confidence level for consistency with the Standard Model in the search 
for f+f" production for different values of m~- and Am. 
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Am 


TO-± (GeV) 

A 1 


{GvV) 


50 


(iO 


70 <S() 


!)() 


94 


signal selection efficiency of tlie general selection at f 89 GeV (%) 


2 


7.9±0.4 


6.6±0.4 


4.9±0.3 


2.3±0.2 


0.5±0.1 


O.liO.O 


3 


21.6±0.7 


21.6±0.7 


20.6±0.6 


19.5±0.6 


15.8±0.6 


13.3±0.5 


4 


34.3±0.8 


33.5±0.7 


33.7±0.7 


33.6±0.7 


31.5±0.7 


32.5±0.7 


5 


41.3±0.8 


43.4±0.8 


42.5±0.8 


41.9±0.8 


43.0±0.8 


44.7±0.8 


10 


61.0±0.8 


64.4±0.8 


63.6±0.8 


63.9±0.8 


63.5±0.8 


65.4±0.8 


20 




75.8±0.7 


77.4±0.7 


78.2±0.7 


79.5±0.6 


78.3±0.7 


(m- 15)/2 




77.1±0.7 


82.2±0.6 


82.7±0.6 


85.8±0.6 


86.6±0.5 


m — 35 


69.7±0.7 


78.3±0.7 


83.2±0.6 


85.9±0.6 


87.4±0.5 


89.0±0.5 


95% CL upper limit on the cross- section times BR'^{xf ~^ 


i^^i)) (fb) 




2 


241.3 


262.3 


330.0 


755.6 


3368.0 


20947.5 


3 


103.8 


89.2 


98.6 


103.8 


88.6 


125.3 


4 


62.3 


60.7 


68.3 


61.5 


80.1 


54.7 


5 


53.0 


44.9 


58.5 


50.3 


63.6 


62.5 


10 


92.2 


82.4 


57.7 


39.0 


39.4 


38.2 


20 




72.1 


64.8 


60.8 


65.3 


40.8 


(m- 15)/2 




81.9 


83.9 


88.2 


61.5 


53.7 


m — 35 


84.6 


87.0 


116.9 


118.5 


70.1 


42.6 


expected upper limit on the cross-section times BIi^{xf — 


i^i>) (fb) 




2 


289.4 


313.9 


387.6 


900.9 


4055.8 


22119.2 


3 


123.1 


117.9 


108.4 


119.4 


134.2 


165.3 


4 


86.7 


82.2 


75.4 


70.9 


77.4 


73.1 


5 


77.5 


68.4 


68.7 


64.7 


59.0 


60.1 


10 


79.5 


64.0 


58.2 


52.1 


47.8 


44.3 


20 




108.5 


84.7 


68.2 


49.7 


38.7 


(m- 15)/2 




121.0 


119.0 


115.5 


108.3 


74.7 


m — 35 


103.8 


130.8 


154.0 


166.4 


143.7 


88.8 


CL for consistency with SM (%) 


2 


63.3 


59.1 


68.3 


70.7 


76.5 


66.1 


3 


61.4 


72.3 


48.1 


56.9 


100.0 


73.5 


4 


80.1 


77.6 


51.8 


53.4 


36.9 


70.2 


5 


82.9 


88.6 


59.0 


72.0 


33.5 


29.5 


10 


24.8 


14.7 


42.1 


76.4 


64.2 


59.2 


20 




80.9 


70.0 


53.1 


15.2 


33.6 


(to - 15)/2 




79.5 


76.0 


69.4 


93.4 


79.0 


m — 35 


63.9 


81.4 


63.6 


72.7 


97.0 


98.9 



Table 6: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 

on the cross-section times BR^{xi — * ^^t'), expected upper limit on the cross-section times 
BR^{Xi (^i^), and the confidence level for consistency with the Standard Model in the search 
ioT xj/^X]^ (2-body decays) production for different values of m^± and Am. The bins without 
entries correspond to values of mp < 35, which are excluded and therefore not considered in 
the analysis. 
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Am 


m-± (GeV) 

Al 


(GeV) 


50 


60 


70 


80 


90 


94 


signal selection efficiency of the 


general selection at 189 GeV (%) 




3 


3.7±0.3 


3.2±0.3 


2.5±0.2 


1.0±0.2 


0.2±0.1 


0.2±0.1 


5 


17.1±0.6 


16.3±0.6 


15.2±0.6 


13.3±0.5 


11.4±0.5 


10.8±0.5 


10 


39.9±0.8 


40.8±0.8 


40.9±0.8 


41.6±0.8 


40.8±0.8 


41.6±0.8 


20 


59.1±0.8 


60.4±0.8 


62.2±0.8 


60.8±0.8 


62.8±0.8 


63.8±0.8 


m/2 


63.3±0.8 


68.7±0.7 


74.3±0.7 


77.9±0.7 


79.3±0.6 


80.3±0.6 


m- 20 


67.3±0.7 


72.6±0.7 


78.0±0.7 


81.2±0.6 


83.9±0.6 


84.0±0.6 


m- 10 


70.1±0.7 


75.9±0.7 


79.4±0.6 


83.0±0.6 


87.0±0.5 


88.1±0.5 


m 


73.2±0.7 


77.4±0.7 


81.2±0.6 


85.2±0.6 


88.1±0.5 


89.5±0.5 


95% CL 


upper limit 


on the cross-section times BFi^{xf ^"^^Xi) (fb) 


3 


355.9 


540.9 


598.9 


1578.0 


5613.2 


11457.9 


5 


111.0 


116.1 


134.6 


114.5 


156.9 


272.9 


10 


74.7 


60.4 


66.5 


61.0 


56.2 


65.5 


20 


91.3 


78.4 


68.0 


64.3 


51.0 


52.8 


m/2 


105.8 


99.7 


82.9 


69.9 


88.0 


121.6 


m- 20 


118.5 


85.2 


91.4 


91.8 


89.2 


89.4 


m- 10 


111.6 


103.7 


104.2 


101.6 


88.1 


60.6 


m 


133.7 


132.4 


117.3 


117.2 


99.6 


57.6 


expected upper limit on the cross- section times BR'^ixf i^uxi) (fb) 


3 


563.5 


668.8 


791.2 


2050.1 


8303.3 


10493.8 


5 


163.4 


154.9 


165.4 


169.6 


203.6 


222.4 


10 


90.4 


82.6 


80.5 


76.1 


76.4 


74.4 


20 


92.0 


81.2 


70.0 


67.3 


63.0 


61.7 


m/2 


100.3 


94.4 


89.9 


81.2 


78.5 


82.5 


m- 20 


103.7 


114.8 


114.9 


113.9 


131.0 


140.5 


m, — 10 


128.2 


130.9 


135.3 


136.9 


166.9 


136.5 


rn 


139.1 


144.8 


irj4.9 


183.3 


179.4 


133.9 


CL for consistency with SM (%) 


3 


94.9 


62.7 


83.3 


76.7 


100.0 


37.5 


5 


82.9 


72.6 


62.5 


88.5 


79.9 


30.5 


10 


62.1 


76.4 


62.6 


68.0 


74.8 


53.0 


20 


41.3 


44.3 


43.7 


46.3 


65.2 


60.2 


m/2 


35.7 


35.7 


49.8 


54.9 


28.6 


8.4 


m- 20 


26.7 


69.0 


65.0 


62.3 


80.5 


84.9 


m- 10 


52.9 


63.4 


66.8 


69.9 


94.2 


99.0 


m 


40.9 


43.9 


64.9 


78.6 


89.8 


99.0 



Tabic 7: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 
on the cross-section times BB?{xf ^^i^Xi), expected upper limit on the cross-section times 
BR^{xf — > l^vXi)-, and the confidence level for consistency with the Standard Model in the 
search for x]''x5"(3-body decays) production for different values of m-± and Am. 



24 



mH+ (GeV) 


45 


55 


65 


75 


85 


94 


signal selection efficiency of the general selection at 189 GeV (%) 


66.6±0.7 


72.3±0.7 


74.6±0.7 


76.9±0.7 


77.6±0.7 


79.7±0.6 


95% CL upper limit on the cross- 


-section times i3i?"^(H* 


^ T^^r) (fb) 


113.6 


84.8 


81.5 


73.2 


77.0 


80.1 


expected upper limit 


on the cross-section times 5i?^(H^ 


-- ^ T^Vr) (fb) 


117.7 


113.2 


114.0 


114.0 


112.3 


111.4 


CL for consistency with SM (%) 


40.8 


70.. S 


7.^).() 


<S3..5 


7(S.7 


7(i.<S 



Table 8: Signal selection efficiency of the general selection at 189 GeV, 95% CL upper limit 
on the cross-section times 5i?^(H='= — > expected upper limit on the cross-section times 

BB?{Y{'^ — > t^Vt), and the confidence level for consistency with the Standard Model in the 
search for H+H~ production for different values of mn-n ■ 
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Figure 1: Distributions of (a) the lepton momentum divided by the beam energy, (b) —qcos6 
and (c) acoUnearity (in radians), for the event sample produced by the general selection at 
y/s — 189 GeV. The data are shown as the points with error bars. The Monte Carlo prediction 
for 4-fermion processes with genuine prompt missing energy and momentum (i~^i> i~T') is shown 
as the open histogram and the background, arising mainly from processes with four charged 
leptons in the final state, is shown as the shaded histogram. In (b) the dashed histogram corre- 
sponds to the distribution expected from smuon pair production, with arbitrary normalisation. 
In (a) and (b) there are two entries per event for events containing two identified leptons. 
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Figure 2: Distributions in the (xmim^max) plane, where x^in x^ax momenta of 

the higher and lower momentum leptou respectively, scaled by the beam energy, for three 
ranges of acolinearity (shown in radians). The distributions are shown for a smuon signal with 
m=45 GeV, Am=45 GeV (left), and for Standard Model Monte Carlo (right). 
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Figure 3: Distributions of acolinearity (in radians) for three example smuon signals. 



28 



OPAL 



-| — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



I 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 4: Distributions of the relative likelihood, L^, for Standard Model Monte Carlo (shaded 
histogram), signal (open histogram) and data (points with error bars), in the analysis for 
smuons with a mass of 80 GeV for a smuon-neutralino mass difference of 60 GeV. 
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Figure 5: Distributions of the background likelihood, Lb, for Standard Model Monte Carlo 
(shaded histogram) and data (points with error bars) for events passing the general selection, 
using all the likelihood variables (a), (b) and (c) show the same information after making 



the initial lepton identification requirements given in Section \4.^ for the selectron and smuon 
searches respectively. In (b), only the variables used in the selectron analysis are used. 
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Figure 6: Contours of the 95% CL upper limits on the selectron pair cross-section times 
BR^ie — i> exi) at 189 GeV based on combining the 183 and 189 GeV data-sets assuming a 
dependence of the cross-section. The kinematically allowed region is indicated by the 
dashed line. The unshaded region at very low Am is experimentally inaccessible in this search. 
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Figure 7: Contours of the 95% CL upper limits on the smuon pair cross-section times 
BE^i^fl fiXi) GeV based on combining the 183 and 189 GeV data-sets assuming 

a dependence of the cross- section. The kinematically allowed region is indicated by the 
dashed line. The unshaded region at very low Am is experimentally inaccessible in this search. 
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Figure 8: Contours of the 95% CL upper limits on the stau pair cross-section times BR^{f — >■ 
rxi) at 189 GeV based on combining the 183 and 189 GeV data-sets assuming a P^/s depen- 
dence of the cross-section. The kinematically allowed region is indicated by the dashed line. 
The unshaded region at very low Am is experimentally inaccessible in this search. 
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Figure 9: Contours of the 95% CL upper limits on the chargino pair cross-section times branch- 
ing ratio squared for xf ^^i^ (2-body decay) at a/s = 189 GeV. The hmits have been cal- 
culated for the case where the three sneutrino generations are mass degenerate. The limit 
is obtained by combining the 183 and 189 GeV data-sets assuming a [3 / s dependence of the 
cross- section. The kinematically allowed region is indicated by the dashed line. The unshaded 
region at very low Am is experimentally inaccessible in this search. 
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Figure 10: Contours of the 95% CL upper limits on the chargino pair cross-section times 
branching ratio squared for xf ^^^^Xi (3-body decay) at y/s = 189 GeV. The hmit is 
obtained by combining the 183 and 189 GeV data-sets assuming a P/s dependence of the cross- 
section. The kinematically allowed region is indicated by the dashed line. The unshaded region 
at very low Am is experimentally inaccessible in this search. 
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Figure 11: The solid line shows the 95% CL upper limit on the charged Higgs pair production 
cross-section times branching ratio squared for the decay — > r^i/r at \/s = 189 GeV. The 
limit is obtained by combining the 183 and 189 GeV data-sets assuming the mH+ and \/s 
dependence oi the cross-section predicted by Pythia. For comparison, the dashed curve shows 
the prediction from Pythia at y/s = 189 GeV assuming a 100% branching ratio for the decay 
— »• t^z/t-. The expected limit calculated from Monte Carlo alone is indicated by the dash- 
dotted line. 
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Figure 12: 95% CL exclusion region for right-handed smuon pair production obtained by com- 
bining the a/s = 183 and 189 GeV data-sets. The hmits are calculated for several values of the 
branching ratio squared for jl^ — >■ /i^Xi ^^^^ indicated in the figure. Otherwise they have 
no supersymmetry model assumptions. The kinematically allowed region is indicated by the 
dashed line. The expected limit for BR^ = 1.0, calculated from Monte Carlo alone, is indicated 
by the dash-dotted line. 
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Figure 13: 95% CL exclusion region for right-handed stau pair production obtained by com- 
bining the = 183 and 189 GeV data-sets. The hniits are calculated for several values of 
the branching ratio squared for t^Xi- The selection efEciency for f^f~ is calculated 

for the case that the decay f" t^X^ produces unpolarised r^. Otherwise the limits have 
no supersymmetry model assumptions. The hatched area shows the region in which the limit 
for BR^ =1.0 can vary if stau mixing occurs (see text). The kinematically allowed region is 
indicated by the dashed line. The expected limit for BB? = 1.0, calculated from Monte Carlo 
alone, is indicated by the dash-dotted line. 
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Figure 14: For two values of tan (3 and /i < —100 GeV, 95% CL exclusion regions for right- 
handed selectron pairs within the MSSM, obtained by combining the y/s = 183 and 189 GeV 
data-sets. The excluded regions are calculated taking into account the predicted branching ratio 
for — >• e^Xi- The gauge unification relation, Mi = ^tan^ 9iyM2, is assumed in calculating 
the MSSM cross-sections and branching ratios. The kinematically allowed region is indicated 
by the dashed line. 
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Figure 15: For two values of tan (3 and fi < —100 GeV, 95% CL exclusion regions for right- 
handed smuon pairs within the MSSM, obtained by combining the \/s = 183 and 189 GeV 
data-sets. The excluded regions are calculated taking into account the predicted branching 
ratio for fl^ — >• /i^Xi- The gauge unification relation, Mi = |tan^ 6*1^^2, is assumed in 
calculating the MSSM branching ratios. The kinematically allowed region is indicated by the 
dashed line. 
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Figure 16: For two values of tan P and < —100 GeV, 95% CL exclusion regions for right- 
handed stau pairs within the MSSM, obtained by combining the \/s = 183 and 189 GeV data- 
sets. The excluded regions are calculated taking into account the predicted branching ratio for 
T^Xi- The gauge unification relation, Mi = |tan^^vF^2, is assumed in calculating the 
MSSM branching ratios. The selection efficiency for f^f~ is calculated for the case that the 
decay f~ —>■ r~Xi produces unpolarised r^. The kinematically allowed region is indicated by 
the dashed line. 
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